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Abstract. We study the diffractive production of high p t photons at HERA. We have 
■ implemented the process as a new hard sub-process in the HERWIG event generator 

^sO , in order to prepare the ground for a future measurement. 

OO ■ 

<n : 

. PACS numbers: 12.38Cy, 13.60.Fz, 13.85.-t 

ON ■ 

5 : 

1. Introduction 

a: 

Ph| One of the cleanest of all diffractive processes is that of diffractive photon production in 
the process ■yp — > 7F where the photon carries a large p t and is well separated in rapidity 
from the hadronic system Y. The process can be measured at HERA [1, 2, 3]. The 
largeness of the transferred momentum — t ~ p\ ^> Aq CD ensures the applicability of 
perturbative QCD. Unlike diffractive meson production this process has the advantage 
that the hard subprocess is completely calculable in perturbation theory. The only 
non-perturbative component resides in the parton density functions of the proton that 
factorize in the usual manner. 

Theoretical interest in this process dates back to the work of [4] where calculations 
were performed in fixed order perturbation theory and to lowest order in <x s . Recent 
work has extended this calculation to sum all leading logarithms in energy, for real 
incoming photons [2] and for real and virtual incoming photons [3]. The cross-section 
for 75 — > 7g can be written 



dpi 167TS 2 ^ ++ ^ ^ ^ 

and we have ignored a small contribution that flips the helicity of the incoming photon. 
The photon-quark CM energy is given by s. To leading logarithmic accuracy [2, 3] 

a ++ = iaa >, E< w Ar rr^r^nw— *w < 2 > 

g y 6 pf J -00 1 + v l [y l + 1/4) 2 -ku 

§ On leave of absence from f . 
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where 

z = — log 4, (3) 
7T pf 

X(v) = 2(*(1) -Re*(l/2 + ii/)) is the BFKL eigenfunction [5], F(z/) = 2(11 + 12z/ 2 ) for 
on-shell photons and there is a sum over the quark charges squared, e 2 . The separation 
in rapidity between the struck parton and the final-state photon is Arj ~ log(s/p 2 ). 

The full photon-proton cross-section is obtained after multiplying by the parton 
density functions: 



da 
dxdp 2 



81 



dpi [) 



H 

and we take the factorization scale \i = p t . 

We have implemented this result in the HERWIG event generator in order to aid 
the experimental measurement of the process. We also note that having done this, 
it is a straightforward procedure to include the high-p t production of vector mesons 
in a similar manner and we intend to do this in the near future. In this paper, we 
compare the HERWIG generated data with theory and discuss the strategy for a future 
measurement. 

In order to speed up the event generation procedure, we used the following 
approximate parameterization: 

G(,) S I- *-*- ^ 2(ll + 12„>^> 

J-oo 1 + v 2 (z/ 2 + 1/4) 2 1XV 
4 52 

« ( ^ + (U)3/2 e 4 - ln2 e(z - 1) + (23.7 + 35^ 3 )0(1 - z) (5) 
which is good to within a few percent over the z-range of interest. 



2. Results 



We show the p t spectrum of the scattered photon in Figure 1. This plot is computed at 
fixed W = 200 GeV (W 2 = s/x) and a fixed a s = 0.2. The choice to fix a s is supported 
by Tevatron and HERA data on gaps between jets and high p t diffractive vector meson 
production [6]. The solid curve shows the theoretical prediction derived directly from 
(4), it is compared to the HERWIG generated data at the parton level (before parton 
showers) and at the hadron level. We used the proton parton density functions of [7] 
(code 155 in PDFLIB [8]). The good agreement between theory and parton level is 
a check that the process is correctly implemented in HERWIG. The systematic shift 
arises because HERWIG ensures that energy and momentum are conserved and that 
final state hadrons are on shell [9]. 

In subsequent figures, we make the typical HERA cuts on the photon energy 
variable, 0.25 < y < 0.75, and on the photon virtuality, Q 2 < 0.01 GeV 2 . Statistical 
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Figure 1. Figure 2. 

The photon p t spectrum. The rapidities of the scattered 

photon and edge of system Y. 



errors are shown corresponding to 44 pb 1 of ep data, typical of that already collected 
by each HERA experiment. We also make a cut yjp < 0.01 where 

(E-p z )i _P^_ a v , g) 
Vf 4- 2£ 7 xW 2 1 ' 

and the sum is over all final state particles excluding the electron and photon. Note 
that yjp can be measured accurately without needing to see the whole of system Y since 
the material lost at low angles does not contribute much to the numerator. As the last 
approximate equality shows, this cut ensures that the rapidity gap between the outoing 
struck parton and the outoing photon is bigger than about 4.5 units (recall that a large 
rapidity gap is a signal of diffractive processes). We have also integrated over all photon 
p t > 2.5 GeV (and used ptmin = 1 GeV in the event generation). 

The effect of varying the yjp cut on the size of rapidity gap can be seen in Figure 2, 
where we plot the rapidity of the scattered photon and the edge of system Y . By not 
requiring to measure system Y in detail it is possible to reach very high rapidity gaps. 
As the plot shows, gaps of 6 to 7 units in rapidity are not uncommon. Note that the 
limited y-range combined with the steeply falling p t spectrum constrain the rapidity of 
the photon to be around r) ~ —2. 

In Figure 3 we show the xjp distribution and in Figure 4 the p t distribution: 

XF 2E p 1 ' 

This variable can be measured to high accuracy. The steep rise at small xjp is driven 
by the BFKL kernel x( u ) i n (2)- I n particular, the dominant conribution comes from 
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Figure 3. 

The xp distribution. 
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Figure 4. 

The pt distribution. 



v « 0, and this leads to 



cixjp pf W 2 \xpj 

where u = (3a s /n)4:\n2 in the LLA. It will be interesting to see to what degree the 
measured xjp distribution follows this power-like behaviour. 



Acknowledgments 



We thank Jon Butterworth and Norman Evanson for their help. 



References 



[1] Ginzburg I F and Ivanov D Yu 1996 Phys. Rev. D54 5523 

[2] Ivanov D Yu and Wusthoff M 1999 Eur. J. Phys. C8 107 

[3] Evanson N G and Forshaw J R 1999 Phys. Rev. D60 034016 

[4] Ginzburg I F, Panfil S L and Serbo V G 1987 Nucl. Phys. B284 685 

[5] Balistky I and Lipatov L N 1978 Sov. J. Nucl. Phys. 28 822 

[6] Cox B E, Forshaw J R and Lonnblad L 1999 JEEP 10 023 

Forshaw J R, 1999 Nucl. Phys. (Proc. Suppl.) 79 311 

[7] Gliick M, Reya E and Vogt A 1995 Z. Phys. C67 433 

[8] Plothow-Besch H, W5051 PDFLIB, 1997.07.02, CERN-PPE 

[9] Cox B E and Forshaw J R 1998 Phys. Lett. B434 133 



